RAOC-T*- 78-158 
maw  Rtport 
July  1978 


CM 


a 

10 


A rUWAHEOTAL  STJOY  OF  THE  ELECTROMASRETIC 
PROPERTIES  OF  ADVANCED  COMPOSITE  MATERIALS 


S ••  * r 3. c ’.»< e Un  tv* r o 1 ty 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  P ACE  (lFh*n  Dmi»  Entered) 

I ([  /R  ^REPORT  DOCUMENTATION  PAGE 


yRADC-yTR-78"15S 
.*■  title  Tina  Sysmul 


S AfTC  READ  INSTRUCTIONS 

BEFORE  COMPLF TING  F O Ril_ 

2 GOVT  ACCESSION  NO.'  » RECIPIENT1*  CATALOG'  NUMBER 


)-k-  Hi-NDAMENTAL  jJJTyDY  OF  THE  ELECTROMAGNETIC 
T^OTERTIES  OF  ^IVANCED ^COMPOSITE  MATERIALS 


I Phase  jtepm/ 

1 Oct  76  - 3 4 Sep  77, 


|7  AuTHORf  ) 


|w.  J.jGajda^^^ 


8.  CONTRACT  OR  GRANT  NUMRERf*) 


|Jf3^2-75-C-^121 


»■  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS  10.  PROGRAM  ELEMENT  PROJECT  T 

AREA  A WORK  UNIT  ^UMBERS 

Syracuse  University  , 61 1 m F /iH  Ih-t 

Syracuse  NY  13210  ^^17^T03  ( « UjJ_  _j 

tl.  CONTROLLING  OFFICE  NAMC  ANO  ADORES5  REPORT  OATE^ 

Rome  Air  Development  Center  (RBCT)  /]  I / 

Griffiss  AFB  NY  13441  I JPPAlTtl 

44 

MONITORING  AGENCY  NAM*  b ADDRESS(7/  dllletent  from  Controlling  Otlicm ) IS.  SECURITY  CLASS,  (of  thin  report) 

Same  UNCLASSIFIED 


15*.  OECL  ASSI  FI  CATION  / DOWN  GRADIN'-. 

schedule 


|'«  Distribution  statement  <o:  mi.  R.ponj 


Approved  for  public  release;  distribution  unlimited. 


lira 


4b#fr#cl  entered  In  Bloch  20,  II  dHtlerent  from  Report) 


SEP  26  1978 


<■  supplementary  NOTES 


R A nr  PlTl  I Qo  t Cf>r>4nnnw>t  !)>••  P r a. a.  a.  - /nnnmt 

v . i\\jj  i . oaaugn  v r\.DL  1 ) 

*Work  done  by  Prof.  W,  Gajda  at  Notre  Dame  University  under  subcontract  with 
Syracuse  University. 


118  KEY  WORDS  fConfin 


old*  II  neceeeery  end  Identify  by  block  number ) 


Advanced  composite  materials 

Permittivity 

Permeability 

Conductivity 


Doping 

Modeling 

Graphite 

Boron 


ASSTR  ACT  (Continue  on  t overee  tide  II  necmeemry  +nd  I dent  1 1)  by  block  number)  — ~~ 

This  report  coveri  an  effort  to  understand,  model  and  modify  the  fundamental 
elt  .trical  properties  of  advanced  composite  materials,  specifically  graphite/ 
epoxy  and  boron/epoxy.  The  electrical  properties  of  graphite  and  boron  fibers 
and  of  epoxy  were  measured.  The  snisotropic  conductivity  of  boron  fibers  was 
studied.  Models  were  developed  which  predicted  the  electrical  properties  of 
Lite  composites — with  reasonable  accuracy  for  some  cases.  The  conductivity  of 
graphite  fibers  was  increased  about  a factor  of  50  (to  lfl^mhos/meter)  by  doping 


do  ,;r7, 1473 


_ UNCLASSIFIED 

EtcuRiTT  clarification  of  thii  Rage  rwh.n  d„,  Enr/u.u; 


A A a H M A 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  Of  THU  RAPt^RTi^i  Dmtm  (nla^j 

wirh  boron;  che  conductivity  of  boron  fiber*  va*  increased  about  a factor  of 
10-*  by  doping  with  carbon. 


i • 


PREFACE 
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Syracuse  University  (Department  of  Electrical  Engineering),  Purdue  University 
(School  of  Electrical  Engineering) » Georgia  Institute  of  Technology  (School 
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(Department  of  Electrical  Engineering)  act  as  prime  contractor  schools  with 
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Air  Force  Academy  (Department  of  Electrical  Engineering),  Air  Force  Institute 
of  Technology  (Department  of  Electrical  Enrineerlng) , and  the  Naval  Post 
Graduate  School  (Department  of  El  cctrical  Engineering)  al«u  p«i  Licipate  in 
the  program. 

The  Post-Doctoral  Program  provides  an  opportunity  for  faculty  at  partic- 
ipating universities  to  spend  up  to  one  year  full  time  on  exploratory  devel- 
opment and  problem-solving  efforts  with  the  post-doctorals  splitting  their 
tine  between  the  cimtomer  location  and  their  educational  institutions.  The 
program  la  totally  customer-funded  with  current  projects  being  undertaken  for 
Rome  Air  Development  Center  (RADC) , Space  and  Missile  Systems  Organization 
(SAMSO),  Aeronautical  System  Division  (ASD) , Electronic  Systems  Division 
(BSD),  Air  Force  Avionlca  Laboratory  (AFAL),  Foreign  Technology  Division  (FTD) , 
Air  Force  Weapons  Laboratory  (AFWL) , Armament  Development  and  Test  Center 
(ADTC) , Air  Force  Communications  Service  (AFCS),  Aerospace  Defense  Command 
(ADC),  Hq  USAF,  Defense  Communications  Agency  (DCA) , Navy,  Army,  Aerospace 
Medical  Division  (AMD),  and  federal  Aviation  Admioliiration  (KAA) . 

Further  information  about  the  RADC  Post-Doctoral  Program  can  be  obtained 
from  Mr,  Jacob  Scherer,  RADC/RBC,  Griffias  AFB,  NY,  13441,  telephone  Autovon 
587-2543,  Coanercial  (315)  330-2543. 
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1)  INTRODUCTION 


Advanced  composite  materials  have  begun  to  be  more  widely  used  as  struc- 
tural and  surface  components  in  modern  aircraft.  Immediate  questions  arise 
concerning  the  effects  of  lightning  upon  composite  structures,  the  shielding 
effectiveness  of  the  materials  against  electromagnetic  interference  and  elec- 
tromagnetic pulse  effects,  the  radar  cross-section  and  precipitation  static 
effects. 

Research  has  shown  that  composites  are  more  likely  than  metals  to  degrade 
structurally  under  simulated  lightning  excitations-'-  and  that  the  electrical 
shielding  effectiveness  of  typical  composites  (graphite/epoxy,  boron/epoxy 
and  Kevlar/epoxy)  is  less  than  that  of  metals  over  the  range  of  frequencies 
DC  to  3 GHz, 2, 3 although  the  shielding  effectiveness  of  graphite/epoxy  ap- 
proaches that  of  metals  for  frequencies  above  100  MHz. 

Partial  solutions  to  the  problems  of  lightning  susceptibility  and  shield- 
ing have  involved  conductive  coatings,  embedded  conductive  grids  and  diverter 
strips . 

However,  a more  fundamental  approach  could  offer  significant  advantages. 
To  anyone  concerned  with  the  problem  of  shielding  the  interior  of  a vehicle 
against  electromagnetic  hazards,  the  possibility  of  a material  with  the  spe- 
cific strength  of  a graphite/epoxy  or  a boron/epoxy  composite  and  the  conduc- 
tivity of  aluminum  is  appealing.  In  order  to  improve  the  shielding  effective- 
ness and  lightning  strike  resistance  of  these  materials,  it  is  natural  to 
consider  ways  in  which  the  electrical  conductivity  and/or  permeability  can  be 
increased  without  significant  decreases  in  specific  strength. 

A necessary  prelude  to  developing  such  materials  is  a determination  of 
the  fundamental  mechanisms  responsible  for  the  electrical  behavior  of  compos- 
ites. Previous  work^  indicates  that  low  frequency  conductivity  depends  on 
fiber  conductivity  and  f iber-to-f iber  contact.  This  leads  to  the  possibility 
of  increasing  fiber  conductivity  by  doping  techniques.  For  example,  poly- 
crystalline boron  is  a semiconductor  and  it  should  be  possible  to  identify 
dopants  which,  when  added  to  the  fibers  in  parts  per  billion  concentrations, 
result  in  orders  of  magnitude  increase  in  fiber  and  composite  conductivity. 
Polycrystalline  carbon  is  less  well  characterized  electrically  but,  in  view 
of  carbon's  position  in  the  fourth  group  of  the  periodic  table,  doping  should 
be  possible  in  this  material  also. 

Permeability  increases  would  improve  shielding  effectiveness  and  could  be 
accomplished  by  doping  composite  fibers  with  atoms  of  high  magnetic  moment. 
Increases  in  mass  density  are  a possible  problem  with  this  approach  and  in- 
creased permeability  would  not  improve,  response  to  lightning. 

Further,  such  doping  of  fibers  should  add  little  to  composite  fabrication 
costs  because  it  could  probably  be  accomplished  by  simply  changing  the  gas 
mixture  in  which  the  fibers  are  formed  (either  by  graphitization  of  a hydro- 
carbon precursor  or  by  chemical  vapor  deposition  of  boron  on  tungsten). 

To  examine  these  ideas,  a research  program  was  initiated  with  the  follow- 
ing three  general  goals: 
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1) 


Characterize  the  electrical  behavior  of  as-manufactured  graphite  and  boron 
fibers  b'  measuring  the  low  frequency  conductivity,  o.  This  t ask  requires 
develop!  ig  reliable  techniques  of  making  ohmic  electrical  contact  to  the 
f ibers. 

2)  Experimentally  and  theoretically  determine  the  fundamental  physics  of  cur- 
rent flow  in  advanced  composite  materials;  in  particular,  if  the  poly- 
crystalline  fibers  are  semiconductor-like,  what  elements  could  be  used  to 
dope  the  fibers  to  change  their  electrical  conductivity? 

3)  Experimentally  investigate  various  doping  methods  to  Increase  the  intrin- 
sic conductivities  of  composites. 

The  following  specific  tasks  were  accomplished  during  the  year.  Tech- 
niques were  developed  for  making  electrical  contact  to  graphite  and  boron 
fibt„..i  and  the  DC  electrical  conductivities  of  fibers  provided  by  three 
manufacturers  (Narmco  and  Hercules  (graphite)  and  Avco  (boron))  were  measured , 
Dit fusions  of  boron  into  graphite  fibers  and  carbon  into  boron  fibers  were 
carried  out  and  the  resultant  changes  in  electrical  conductivity  were  measured. 
The  high  field  non-linear  thresholds  of  the  fibers  were  measured  in  order  to 
assess  failure  mechanisms  under  high-  level  electrical  excitation.  Models  were 
formulated  to  allow  the  prediction  of  composite  conductivities  based  on  the 
conductivities  of  the  constituent  fibers  and  matrix.  It  was  theoretically 
concluded  that  attempts  to  increase  the  magnetic  permeabilities  of  advanced 
composite  materials  would  result  in  unacceptable  decreases  in  specific 
strength  and  this  approach  was  not  investigated  experimentally.  Such  a spe- 
cific strength  decrease  will  result  from  the  mass  density  increase  as  rela- 
tively dense  high  permeability  materials  are  added.  For  example,  loading  the 
composite  with  10  volume  percent  high  permeability  iron-nickel  alloy  would 
increase  the  composite  density  by  402  - a change  which  would  be  unacceptable 
in  many  applications  although  the  relative  permeability  would  be  20. 

The  next  section  of  this  report  develops  simple  models  of  current  flow  in 
advanced  composites.  The  large  differences  in  the  conductivities  of  graphite/ 
epoxy  and  boron/epoxy  are  demonstr;  ted  and  the  central  role  played  by  the 
fibers  in  determining  overall  material  conductivity  is  illustrated.  These  re- 
sults indicate  that  increasing  fiber  conductivity  will  indeed  lead  to  high 
composite  conductivities.  Although  such  materials  would  have  to  be  thoroughly 
tested,  a comparison  of  the  responses  of  metals  and  presently  manufactured 
advanced  composites  to  electromagnetic  hazards  indicates  that  enhanced  com- 
posite conductivities  would  lead  to  reduced  electromagnetic  vulnerability. 

Section  3 introduces  the  theory  of  impurity  diffusion  in  cylindrical 
geometiles  germane  to  fibers. 

The  experimental  investigations  of  as-manufactured  graphite  and  boron 
fibers  are  discussed  in  Sections  4 and  5 respectively.  Techniques  used  to 
characterize,  clean  and  form  electrical  contacts  are  described.  The  DC 
electrical  conductivities  and  anisotiopies  associated  with  the  fibers  are 
determined . 

The  experimental  techniques  used  to  introduce  impurities  into  the  fibers 
are  detailed  in  Section  6.  Thermal  diffusion  was  used  for  both  boron  and 
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graphite  while  the  intercalation  of  super-acid  radicals  Into  graphite  was 
also  investigated.  The  results  of  these  processes  a'e  given  in  Section  7. 

The  last  section  summarizes  the  progress  made  and  indicates  the  direc- 
tions along  which  the  research  will  continue. 

2)  CURRENT  IN  COMPOSITES 

The  materials  considered  in  this  report  are  inhomogeneous,  consisting  of 
arrays  of  fibers  (each  fiber,  in  general,  electrically  anisotropic)  in  a 
matrix  (electrically  isotropic  for  all  materials  presently  under  investiga- 
tion). The  theoretical  situation  is  worsened  by  the  lack  of  periodicity  which 
precludes  the  use  of  translational  invariance  symmetry  as  commonly  applied  in 
ths  theory  of  crystalline  solids. 

It  is  possible  to  identify  two  major  lines  of  analysis  of  such  mater ial 5 ■ 
The  first  may  be  termed  the  statistical  study  of  inhomogeneous  materials  and 
has  been  developed  primarily  by  individuals  studying  problems  of  heat  flow  in 
disordered  media  (thermal  insulation,  for  example).  The  basic  theme  is  the 
determination  of  upper  and  lower  bounds  of  the  transport  coefficients  (thermal 
conductivity,  electrical  conductivity,...),  in  terms  of  the  properties  of  the 
two  phases  and  the  volume  fractions  of  each  phase.  Although  such  models  have 
been  developed  primarily  for  solids  consisting  of  spheres  dispersed  in  a ma- 
trix, extensions'  have  been  published  for  dispersions  of  fibers.  It  must  be 
noted  that  the  theoretical  analyses  do  not  assume  any  touching  between  the 
fibers  but  rather  assert  that  euch  fiber  is  surrounded  by  the  matrix  material. 

We  begin  by  examining  the  applicability  of  these  bounds  to  the  case  of 
graphite  fibers  in  either  an  epoxy  or  a polyimide  matrix.  It  is  most  con- 
venient to  examine  the  ratio  of  the  upper  bound  (o^)  to  the  lower  bound  (o^). 
The  general  expression  takes  the  form  5 

°U/0L  “ a~ljl  + af  "AJ  / If  + a (1-f)  + B] 
with  A - . 5af  (1-f)  / [f  + 20^1-f)2  - 2G2f2] 

B - .25  of(l-f)  / (f2G2  - (1-f) 2Gjl ] 

where  a Is  the  ratio  of  fiber  to  matrix  conductivities,  f the  volume  frac- 
tion of  the  fibers  and  G^,  G2  are  empirical  shape  factors  (see  below). 

In  writing  the  above  equation,  available  experimental  evidence  indicat ing 
that  fiber  conductivities  are  several  orders  of  magnitude  greater  than  matrix 
conductivities  was  Invoked.  As  long  as  neither  f nor  1-f  closely  approaches 
aero  (a  condition  which  is  satisfied  in  both  typical  graphite/epoxy  and  graph- 
ite/polylmide  samples),  the  bounds  ratio  (0  /0  ) simplifies  further  and  numer- 
ical values  for  the  shape  factors  may  also  Be  Substituted, 

G1  " G2  " 0-25 
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For  a sain] lc  of  graphite/epoxy  or  graphite/poly iroide  representative  of 
the  materials  used  in  aircraft,  the  following  approximate  numerical  values 
apply. 

6 ”8 

f » 0-7,  Cj  » 10  mhos/meter,  “ 10  mhos/meter 

The  bounds  may  bs  evaluated  directly  as 


— - 4.41(1010) 

\ 


This  value  is  significantly  greater  than  the  experimentally  measured 
values , ^ 


It  may  be  concluded  that  these  Beran-Si lnutzer  bounds  are  not  suf f ic ien t - 
ly  tight  to  realistically  predict  the  electrical  properties  of  graphite/resin 
composites.  Although  the  bounds  are  valid,  thev  are  more  than  ten  orders  of 
magnitude  apart  and  are  of  little  value  in  predicting  the  conductivities  to 
be  expected  in  a given  composite.  These  results  imply  that  the  basic  diver- 
gence between  this  theory  and  experiment  involves  the  theoretical  neglect  of 
f iber-to-f iber  contact.  This  contact  is  a dominant  feature  in  configurations 
in  which  current  is  directed  orthogonal  to  the  fiber  axis  and  cannot  be  ig- 
nored if  useful  bounds  are  to  be  derived  theoretically. 

The  second  approach  to  the  analysis  of  the  electrical  properties  of  com- 
posites rests  upon  deterministic  models  in  which  one  attempts  to  solve  for 
the  precise  distributions  of  currents  in  a sample  excited  by  an  external  volt- 
age and,  in  this  way,  determine  an  effective  conductivity  for  the  material. 
Such  models  become  intractable  when  large  numbers  of  fibers  must  be  considered 
but  insight  can  be  obtained  into  the  basic  electrical  behavior  by  examining  a 
few  simple  models. 

We  begin  our  development  of  such  models  of  electrical  conductivity  in 
advanced  composites  by  considering  the  most  simple  system — the  unidirectional, 
aingle-ply  system.  If  care  if  taken  to  abrade  any  epoxy-rich  regions  from 
the  upper  and  lower  surfaces  of  such  a sample,  we  find  that  the  conductivity 
may  be  characterized  by  two  quantities:  o^,  the  value  associated  with  current 
along  the  fiber  axis  and  o , the  value  associated  with  current  orthogonal  to 
the  fiber  axis  (either  from  top  to  bottom  or  in  the  sample  plane  at  right 
angles  to  the  fibers). 

We  use  the  symbol  f to  denote  the  volume  fraction  of  fibers  in  the 
unidirectional  sample,  o for  the  fiber  conductivity  and  o for  the  epoxy 
matrix  conductivity.  Fof  current  in  the  longitudinal  direction  (along  the 
fibers),  the  total  sample  conductance  (Figure  1)  is  simply  the  summation  ol 
the  conductances  of  the  individual  fibers  plus  the  conductance  of  the  epoxy. 
This  may  be  written  as 
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MODEL  KOH  DETEKMI NATION  OF 

FIGURE  i 

Although  minor  deviations  from  this  result  are  observed  in  some  unidirec- 
tional Simples  primarily  as  a result  oi  breaks  in  fibers  and  skewing  of  fibers 
from  the  longitudinal  axis,  it  yields  acceptable  results.  For  example,  the 
DC  conductivity  of  'typical'  graphite  fibers  has  been  measured''  as  2(10/',)mhos/m. 

(Section  4).  For  unidirectional  samples  with  a fiber  volume  fraction  of  0.6, 
the  above  equation  predicts  the  longitudinal  conductivity  to  be 

oL  - 2 ( 1 04 ) (0. C)  + o (0.4) 

m 1.2(10^)  mhos/m, 

The  second  term  is,  of  course,  negligible  when  the  matrix  is  epoxy.  This 
value  is  in  excellent  agreement  with  low  frequency  experimental  measurements  1 

{DC  to  S MH?)>  which  have  been  in  the  range  cf  1-2(10  } mhoa/m  for  uuiuiiec-  | 

tional  samples.  I 


1 
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Th°  calculation  of  the  transverse  conductivity  oT  1 r.  significantly  more 
difficult.  In  order  to  develop  insight  into  this  problem,  consider  a simple 
composite  composed  of  two  fibers  with  rectangular  cros9-sections,  which  touch 
each  other  via  a rectangular  "neck."  This  model  is  shown  in  Figure  2, 


DIRECTION  OF  CURRENT 
— * 


MODEL  FOR  CALCULATION  OF  oT 


FIGURE  2 


Consider  a slab  thick,  depth  W and  length  L)  of  this  sample  as 
sketched  In  Figure  2 . (u  conductance  of  this  section  can  be  calculated 

directly  as 

AG  = o^WAx/L  for  - y f 3 i s < j - P 
AG  = o^^WAx/to^a  + (^(L-oO)  otherwise 


The  various  geometric  terms  are  defined  in  Figure  2.  a,  6 are  used  to 
define  the  'neck'  connecting  the  two  fibers  as  shown  in  the  figure. 

The  total  conductance  is 


G(+t/2) 
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dG  - °lWt  + 
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_ 
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Th' - result  may  be  rewritten  in  terms  of  a transverse  conductivity  o^, 
defined  as 


G *=  a^tW/L, 
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■■’sm****^,  ■ ■*.-»«* i'1  "4'  


which  implies  that 


a.Oj(2S/t) 

°T  “ °1  (1  “ 26/0  + 0l(a/L)  + o2(l  - a/L) 

When  3 ■+  t/2  and  a L,  the  model  simplifies  to  pure  matrix  and  the  con- 
ductivity takes  the  form 

0T  * °i  U-  ~ 1)  + (°i  + °2  " °2‘ 

When  $ ■+  0 and  a 0,  the  model  simplifies  to  pure  fiber  ana  the  conduc- 
tivity is 


a 


T 


as  expected. 


The  general  expression  for  a is  of  interest  because  of  the  information 
it  contains  about  the  situation  in  which  the  fiber'  io  not  touch.  This  is 
described  by  the  geometric  condition  3 - t/2,  and  the  equation  for  trans- 
verse conductivity  reduces  to 

CT  " . 


The  low  probability  of  non-touching 
is  demonstrated  by  experimental  results, 
graphite/epoxy  has  been  measured  to  have 
quantity  L is  approximately  2(10“^)m. 
epoxy  matrix  to  be  10_®  mhos/m  yields 


fibers  in  graphite/epoxy  materials 
The  transverse  conductivity  for 
a value  of  200  mhos/m. , while  the 
Estimating  the  conductivity  of  the 


1(10-15)  m. 


Such  a small  separation  is,  for  all  meaningful  purposer,  equivalent  to 
contact  between  fibers  and  it  may  be  concluded  that,  in  a transverse  sample 
of  graphite/epoxv , f iber-to-f iber  contact  is  essential  to  describe  the  conduc- 
tion of  charge  in  directions  orthogonal  to  the  fiber  axis. 


It  is  important  to  note  that  such  fiber  contact  does  not  occur  in  boron/ 
epoxy  and,  in  these  materials,  the  flow  of  charge  across  fibers  is  limited  by 
the  conductivity  of  the  epoxy  matrix.  In  unidirectional  samples  of  boron/ 
epoxy  ^ , a ” 30  rahoB/m  and  a ■ 2(10~°)  mhos/m. 


In  composites  which  have  an  insulating  matrix,  two  factors  determine  the 
conductivities  of  unidirectional  samples;  the  fiber  conductivity  and  the 
degree  of  fiber  contact.  For  multiple-ply  lay-ups,  a simple  model^  describ- 
ing the  overall  conductivity  o is  given  by  the  equation 


o 


N 


N 
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where  N is  the  number  of  layers  and  o is  the  conductivity  in  the  direc- 
tion of  interest  of  the  fth  layer.  This1equation  follows  from  the  assumptions 
that  the  conductances  of  the  layers  may  be  treated  as  effectively  in  parallel 
and  that  each  ply  has  the  same  thickness. 

All  of  these  results  pertain  to  situations  in  which  current  is  moving 
between  opposites  faces  of  rectangular  slabs  of  the  material.  Physical  sit- 
uations of  interest  often  display,  of  course,  more  complex  geometries.  In  the 
case  of  lightning  attachments  or  precipitation  static,  the  current  is  injected 
into  a series  of  points  along  the  surface  and  then  distributes  itself  (via  a 
relaxation  process)  through  the  three-dimensional  aircraft  fuselage.  Although 
detailed  modeling  of  the  process  is  difficult,  it  is  clear  that  the  above 
description  will  hold  qualitatively.  There  will  be  differences  in  various 
formulas  as  a result  of  the  geometric  complexity  but  the  basic  conclusion  that 
fiber  contact  dominates  the  observed  composite  material  conductivity  will  hold. 

The  introduction  of  advanced  composite  materials  causes  a large  number 
of  problems  for  individuals  concerned  with  the  electromagnetic  compatibility 
and  survivability  of  defense  systems.  A significant  effort  has  already  been 
expended  to  predict  electromagnetic  scattering  and  penetration  for  all  metal- 
lic structures.  The  applicability  of  this  research  to  advanced  composite 
materials  is  not  known.  If  fundamental  physical  changes  in  composites  could 
be  made  so  that  their  electrical  properties  more  closely  approach  those  of 
metals,  existing  scattering  and  penetration  algorithms  could  be  applied  to 
composites  with  improved  accuracy. 

The  models  of  this  section  indicate  that  this  goal  could  be  approached 
in  one  of  two  ways:  either  increase  the  conductivity  of  the  fibers  or  in- 
crease the  fiber  contact.  The  second  path  is  ruled  out  by  mechanical  require- 
ments. The  resistance  of  composites  to  microbuckling  is  primarily  determined 
by  the  fiber  contact.  As  fiber  contact  increases,  resistance  to  microbuckling 
decreases  and,  in  general,  these  materials  are  fabricated  so  as  to  minimize 
fiber  contact.  Attempts  to  increase  conductivities  by  increasing  fiber  con- 
tact would  lead  to  an  unacceptable  degradation  in  mechanical  behavior. 

This  forces  us  to  examine  fiber  conductivities  more  carefully  and  develop 
techniques  for  increasing  the  conductivities  without  an  accompanying  decrease 
in  specific  strengths.  This  implies  that  the  basic  mechanical  strengths  of 
the  fibers  must  be  unchanged  and,  at  the  same  time,  their  densities  should  not 
increase.  These  requirements,  in  conjunction  with  processing  difficulties, 
rule  out  plating  the  fibers  with  metallic  sheaths. 

The  model  used  above  to  calculate  is  the  simplest  possible  in  that  it 
involves  two  rectangular  fibers  with  a rectangular  neck  between  them.  It  is 
of  interest  to  examine  a more  realistic  model — che  random  fiber  model.  The 
transverse  conductivity  is  to  be  calculated  for  current  through  a series  of 
parallel  planes.  The  rectangular  fibers  are  distributed  across  each  plane 
with  a probability  f,  the  macroscopic  volume  fraction  of  the  fibers.  For  a 
sample  composed  of  two  such  planes  (Figure  3),  the  transverse  conductivity  is 
calculated  under  the  assumption  that  the  boundaries  of  the  two  layerc  arc 
equipotentials. 


i 
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RANDOM  FIBER  MODEL 
FIGURE  3 

At  any  lateral  point  X,  there  are  four  possibilities  concerning  the 
material  in  the  z direction: 

2 

i)  both  layers  are  fiber;  probability  f , 

ii)  the  first  layer  is  fiber,  the  second  matrix;  probability  f(l-f) 

iii)  the  first  layer  is  matrix,  the  second  fiber;  probability  (l-f)f, 

2 

iv)  both  layers  are  matrix;  probability  (1-f)  . 

The  total  current  flow  is  the  sum  of  the  currents  associated  with  each 
of  the  four  cases.  For  case  x, 

Jf  *=  = o^Aifj/t) 


because  the  electric  field  is  uniform  through  the  two  layers  and  has  a value 
Aif>/t  where  A $ is  the  total  potential  drop  and  t/2  is  the  thickness  of  a single 
layer.  (Figure  3.  ) 


In 


type  ii  situations,  the  electric  field  in  the  fiber,  E^ , and  in  the 


matrix,  E^,  are  determined  from  the  conservation  of  current  density. 


'1  °1^1  “2 

and  the  recognition  that  the  total  potential  drop  across  the  two  layers  isAtf>, 
E1(t/2)  + (t/2)  = A* 


°2E2 


J 

5 

j 

* 

i 
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Solving  these  latter  two  equations  simultaneously  yields 


Similarly,  for  the  other  two  cases. 


Jin  ' 2 


Jiv  ” °2 


The  total  current  density  is  found  by  weighting  each  of  the  four  individual 
probabilities  and  summing 

J - + f(l-£)J11  ♦ U-£)£JIU  + U-£)2Jiv 

Recognizing  that  J - o^E  and  E » A$/t  yields 

°T  " °lf2  + 4(V^)f(1_f)  + 02(1-f)2 

For  >>  a 2 (the  situation  with  graphite/resin  composites),  this  re- 
duces to 

oT  ■ a^f2  + 4o2f(l-f)  + o2(l-f)2 

As  long  as  f doesn't  approach  zero,  it  may  be  further  concluded  that 


In  graphite/epoxy  samples  with  f ■ 0.6,  direct  measurement  indicates  that 

4 

oT  - 200  mhos/meter  and  ■ 2(10  ) mhos/meter 
while  the  above  equation  predicts  that  - 7200  mhos/meter. 
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If  the  basic  random  plane  model  is  extended  from  two  to  three  layers, 
the  method  of  analysis  remains  the  same  but  the  algebraic  complexity  in- 
creases. Again,  for  >>  a and  f not  near  zero,  the  transverse  conduc- 
tivity may  be  written  as 


Extending  this  reasoning  produces  the  following  table 


Number  of 

Layers  in  Model  CT^°1  CT  (m'10s/ra,) 

1 f 1.2(104) 

2 £2  7.2(103) 

3 f3  4. 3(103) 


9 f9  2.0(102) 

10  f10  1.2(10") 

11  f11  7.3(10) 


4 

The  above  values  for  were  calculated  using  f=0.6  and  o^  *=  2(10  ) 

mhos/m. 


For  a single  ply  of  graphite/epoxy,  the  nominal  thickness  is  0.005  inches 
while  the  fibers  have  diameters  of  0.0003  inches.  This  implies  that  the  ran- 
dom plane  model  should  contain  approximately  17  layers  (.005/. 0003  = 16.67). 
However,  the  leading  term  converges  to  the  experimental  value  for  only  ten 
layers.  However,  with  f-O.TS,  the  random  plane  model  yields  a value  of 
close  to  th<  experimental  value  with  the  number  of  planes  equal  to  17. 

This  strong  dependence  upon  the  number  of  planes  in  the  model  results 
from  the  use  of  only  the  leading  term  in  the  full  expression  for  o^.  As  the 
number  of  planes  increases,  the  calculation  of  terms  other  than  the  leading 
term  becomes  quite  complex  because  one  has  to  count  all  possible  fiber-to- 
fiber  paths  through  the  sample, 

A FORTRAN  program  has  beer,  written  to  count  such  paths  but  it  io  not,  at 
present,  operating  correctly.  The  algorithm  does  not  find  all  paths  and 
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results  in  incorrect  values  of  conductivity.  The  program  details  are  not 
included  in  this  report  and  this  aspect  of  the  research  is  continuing. 

Although  the  models  examined  in  this  section  have  been  simple,  they  lead 
to  an  important  conclusion.  The  intrinsic  conductivity  of  graphite/resin  and 
boron/resin  advanced  composite  materials  can  be  increased  by  increasing  the 
fiber  conductivities. 


3)  ELEMENTS  OF  MICROSCOPIC  FIBER  CONDUCTIVITY  AND  IMPURITY  DIFFUSION 

Given  the  desirable  goal  of  increasing  fiber  conductivity,  It  is  neces- 
sary to  observe  two  implicit  constraints.  If  the  enhanced  conductivity  fi- 
bers are  to  be  useful  in  advanced  composite  materials,  it  is  necessary  that 
any  process  used  to  increase  fiber  conductivity  must  not  decrease  the  specific 
strengths  of  the  fibers  to  any  significant  degree. 

To  gain  insight  into  possible  methods  for  attaining  these  goals,  the 
equations  for  the  conductivity  a of  a semiconductor  provide  a starting  point. 


o = q(p^n  + u^p) 

-19 

where  q = 1.6(10  ) coulombs, 

u , u,  arc  the  electron,  hole  mobilities  and 
e n 

n,  p are  the  densities  of  "free"  electrons, 
holes  per  unit  volume. 

In  the  case  of  a metal,  the  hole  concentration  is,  of  course,  zero. 

The  electron  and  hole  mobilities  essentially  measure  the  ease  with  which 
charged  particles  move  through  the  material  under  the  influence  of  an  exter- 
nal electric  field.  The  mobilities  depend  primarily  upon  the  atomic  species 
of  which  the  solid  is  composed  as  well  as  the  crystalline  perfection.  In 
osrt9r3l y 3.  hi^hsr  rss  of  cr,r£tsl ii. ns  ordsr  in>nli6s  3l  I(?;r°6r  vsluc  of 
mobility. 

These  remarks  lead  to  the  observation  that  the  conductivities  of  the 
graphite  and  boron  fibers  used  in  advanced  composites  could  be  increased  by 
more  careful  attention  to  the  processes  whereby  the  fibers  are  formed 
(graphitization  of  a hydrocarbon  and  chemical  deposition  of  boron  on  tung- 
sten) to  insure  a higher  degree  of  crystalline  order.  Such  efforts  have  not 
been,  to  our  knowledge,  made. 

Increases  in  electron  and  hole  densities  can  result  In  large  increases 
in  conductivities  and  orders  of  magnitude  increases  result  from  relatively 
small  additions  of  the  proper  impurity  atoms.  In  a typical  semiconductor, 
an  addition  of  0.1  Darts  ner  million  of  a donor  imouritv  increases  the  con- 

. . -I  -i 

ductivity  by  over  five  orders  of  magnitude.  . 
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Although  future  confirming  experiments  will  be  necessary,  it  seems  likely 
that  such  a small  impurity  addition  to  the  fibers  will  not  have  a major  effect 
upon  either  the  mechanical  strength  or  the  mass  density. 

This  process  of  adding  small  amounts  of  electrically  active  impurity 
atoms  (either  donors  or  acceptors)  to  semiconductors  is  known  as  "doping." 

It  is  widely  used  in  fabricating  solid  state  devices  in  silicon,  germanium 
and  compound  semiconductors. 


The  next  germane  question  is  "Are  graphite  and  boron  fibers  semiconduc- 
tors?". Available  evidence  indicates  that  both  materials  are  semiconductors 
in  their  normal  composite  forms.  Semiconductor  behavior  has  been  observed 
in  polycrystalline  specimens  of  both  graphite^  and  boron  ■ 


One  caveat  is  essential  at  this  point.  As  impurity  atoms  are  introduced 
into  a host  solid,  their  degree  of  electrical  activity  depends  upon  the  crys- 
talline perfection  of  the  host  medium.  Given  the  high  degree  of  atomic  dis- 
order present  in  the  fibers  used  in  advanced  composite  materials,  it  is  highly 
probable  that  the  doping  process  will  not  be  as  efficient  as  it  is  in  single 
crystal  silicon.  This  inefficiency  will  lead  to  limits  upon  the  amounts  of 
conductivity  enhancement  attainable  in  both  graphite  and  boron  fibers. 


The  most  natural  method  of  adding  impurities  to  fibers  utilizes  the  ther- 
mal diffusion  of  impurity  atoms  under  the  influence  of  a concentration  gradi- 
ent. Before  describing  the  experimental  aspects  of  Impurity  diffusion,  the 
basic  elements  of  linear  diffusion  theory,  will  be  reviewed. 

The  first  order  theory  assumes  that  the  impurity  atom  flux  J (atoms/cm  - 
sec)  is  linearly  related  to  the  spatial  gradient  of  the  impurity  concentration 
N (atoms/cm^).  The  relationship  between  J and  N is  called  Fick's  First  Law10 
and  is  written 


J = -DVN 


where  D is  the  diffusion  coefficient  (cm^/sec)  and  V the  spatial  gradient 
operator.  The  diffusion  coefficient's  numerical  value  is  usually  obtained 
from  experiment. 

Invoking  the  conservation  of  impurity  atoms  results  in  Fick's  Second 
Lawll  which  may  be  written  as 


~ = V • DVN 

0 L 

where  t is  the  time. 

The  problem  of  interest  involves  the  diffusion  of  impurity  atoms  into 
fibers.  For  analytical  ease,  the  fibers  are  modeled  as  being  infinitely  long 
so  that  end  effects  may  be  neglected  and  the  fibers  are  also  assumed  to  have 
circular  cross-sactional  areas.  Given  that  the  source  of  impurity  atoms  is 
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uniformly  distributed  either  on  or  near  the  fiber  surface,  it  follows  that 
the  distribution  of  the  diffused  impurity  species  (N(r,t))  will  be  independ- 
ent of  the  polar  angle.  The  impurity  concentrat ion  depends  only  upon  the 
radial  distance,  r and  the  time,  t.  In  this  case,  Fick's  Second  Law  takes 
the  form 


3N  1 3_  9N 
3t  “ r 3r  rU  ar‘ 


If  the  initial  distribution  of  impurities  in  the  fiber  is  known  and  the 
source  used  to  add  impurities  to  the  fibers  during  diffusion  is  suitably 
characterized,  solutions  to  this  partial  differential  equation  will  predict 
the  resultant  impurity  profile  for  a diffusion  at  a given  temperature  T and 
time  t.  The  coefficient  D is  a measure  of  the  relative  ease  with  which 
impurity  atoms  move  through  the  host  material.  The  larger  the  value  of  D, 
the  faster  the  impurity  atoms  diffuse.  Although  diffusion  theoretically  oc- 
curs in  any  physical  system  in  which  a non-zero  spatial  gradient  of  impurity 
concentration  exists,  the  diffusion  mechanism  is  so  slow  at  room  temperature 
in  solids  as  to  be  negligible.  As  a consequence,  any  intentional  diffusion 
processing  is  carried  out  at  elevated  temperatures.  A suitable  range  of 
temperatures  for  a given  system  may  be  predicted  from  a knowledge  of  the  dif- 
fusion coefficient  D. 

A solution  to  Fick's  Second  Law  has  been  found11  in  closed-form  for  the 
situation  in  which  a cylinder  of  radius  P.  is  initially  impurity  free  and 
exposed  to  a constant  surface  impurity  concentration  during  an  elevated  tem- 
perature diffusion  cycle  for  a time  t . It  is  assumed  that  the  diffusion 
coefficient  D is  a constant.  The  solution  is 

? . - 2t 

N(r,t)  = N (1  - - l e °n  (J  (ra  ) / J 'Ra  ))] 
o K - o n -L  n 

n=l 


where  J 

o 


is  the  Bessel  function  of  the  first  kind  of  order  zero, 
the  Bessel  function  of  the  first  order  and 


a are  the  roots  of  J 
n o 

Numerical  evaluation11  of 

N(0,to)  * N (R.t^)  - 

D t - 
o 


(Ra  ) - 0. 
n 

this  formula  reveals  that 

N to  within  1%  for 
o 


Although  little  information  could  be  found  regarding  the  diffusion  co- 
efficients of  impurities  of  interest  in  graphite  and  boron,  estimates  of  the 
time  required  to  produce  a uniform  impurity  distribution  may  be  made  using 
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values  o 
D - lO'3 
time  t 

o 


f D for  electrically  active  Impurities  In  silicon.  For  example, 
* cm^/see,  for  phosphorous  diffusing  in  silicon  at  1200°  C.  The 
required  to  produce  a nearly  uniform  impurity  concentration  is 


R'VD 


12  2 
10  R 


For  graphite  fibers, 

i 

R = 3(10  cm 

A 

and  t - 9(10  ) sec  * 25  hours 

o 

This  time  may  be  reduced  significantly  by  increasing  the  temperature  at 
which  the  impurities  are  diffused  since  D has  an  exponential  dependence  on 
the  temperature  T,  i.e.. 


D(T)  » Dq  exp  (-Q/T)  (D^,  Q constants) 

In  summary,  thermal  diffusion  is  a pr.ocess  whereby  impurity  atoms  move 
into  a solid  under  the  influence  of  an  impurity  concentration  gradient.  There 
is.  on  the  average,  a net  flux  of  impurity  atoms  from  regions  of  high  concen- 
tration to  regions  of  low  concentration.  This  aggregate  motion  is  a conse- 
quence of  the  random  walk  motion  of  the  impurity  atoms.  This  process  is  well 
characterized  theoretically  and  has  been  widely  used  experimentally  to  intro- 
duce impurities  into  semiconductors.  For  these  reasons,  high  temperature 
diffusion  has  been  chosen  as  the  method  to  add  selected  impurities  to  graphite 
and  boron  fibers. 

It  is  evident  from  the  discussion  of  this  section  that  little  quantita- 
tive information  ie,  at  present,  available  concerning  impurity  diffusion  in 
graphite  and  boron.  Host  importantly,  numerical  values  for  the  diffusion 
coefficients  as  a function  of  temperature  for  various  impurities  in  graphite 
and  boron  ate  not  known.  The  experimental  work  described  in  the  following 
sections  contains  no  quantitative  analysis  of  the  impurity  distributions  re- 
sulting from  the  diffusion  cycles.  The  experiments  were  intended  to  assess, 
in  a preliminary  way,  the  feasibility  of  modifying  fiber  conductivities  via 
impurity  diffusion.  In  follow-on  work,  secondary  Ion  mass  spectroscopy  will 
be  used  to  obtain  detailed  impurity  profiles  resulting  from  specific  temper- 
ature-time diffusion  cycles  and,  from  this  information,  specific  numerical 
values  for  diffusion  coefficients  will  be  derived.  Such  detailed  information 
was  not  a goal  of  this  year's  research  activities. 

In  the  next  sections,  the  physical  and  electrical  properties  of  as-grown 
fibers  are  examined.  In  addition,  techniques  used  to  make  repeatable,  elec- 
trical contacts  to  both  graphite  and  boron  fibers  are  discussed  in  detail. 
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4)  GRAPHITE  FIBERS 


There  has  been  a significant  amount  of  research  conducted  concerning  the 
response  of  graphite/epoxy  composites  to  laboratory  simulations  of  electro- 
magnetic excitations  ’ ’ Such  material  behavior  must  be  thoroughly  under- 

stood before  aircraft  made  primarily  of  composite  materials  can  be  safely  de- 
signed and  flown.  The  research  to  date  has  been  empirical  and  has  resulted 
in  the  conclusion  that  protective  metal  thin  films  or  mesh  are  necessary  to 
prevent  severe  structural  degradation  from  lightning  attachment.  Some  con- 
trary opinions1^  have  also  been  expressed. 

A complete  model  of  graphite/epoxy  response  to  electromagnetic  excita- 
tion invol\ es  several  complex  elements.  Pre-requi3ltes  include  a knowledge 
of  the  low-level  (linear)  response,  current-voltage  characteristics  of  the 
graphite  fibers,  the  high-level  (non-linear)  response,  a statistical  knowledge 
of  fiber  arrangements  and  some  description  of  f iber-to-f iber  contact,  solu- 
tions detailing  the  distribution  of  currents  in  fiber/epoxv  structures  and  a 
thermal  transient  analysis  of  such  systems.  The  first  task  is  to  examine  the 
electrical  responses  of  the  fibers.  The  low  frequency  conductivity  of  uni- 
directional graphite/epoxy  composites  has  been  reported^1  and  the  need  for 
careful  sample  preparation  and  standardizat  ion  has  been  discussed.  “ In  this 
section,  we  examine  the  linear  current-voltage  characteristics  of  graphite 
fibers  as  well  as  the  threshold  electric  fields  and  current  densities  asso- 
ciated with  failure. 

Preliminary  investigation  revealed  that  fibers  produced  by  various  com- 
panies had  different  conductivities.  To  simplify  the  following  discussions, 
the  work  reported  involved  only  one  fiber  type — Thorne!  T300  as  used  in 
Narmco  5208  pre-preg  tapes. 

Individual  graphite  fibers  have  radii  on  the  order  of  a few  micrometers 
and  are  difficult  to  handle.  In  general,  it  was  convenient  to  immerse  a seg- 
ment of  graphite  tow  in  methanol  and  use  teflon  coated  tweezers  to  separate 
Individual  fibers  from  the  bundle.  Microscopic  observation  and  a steady  hand 
during  this  procedure  improved  the  success  rate.  The  fibers,  in  lengths  from 
two  to  four  centimeters,  were  mounted  on  pre-cleaned  glass  microscope  slides, 
mis  was  accump-Lisutru  oy  applying  eitctncany  conductive  silver  paint  to  tne 
ends  of  the  fibers.  When  dry,  this  paint  holds  the  fiber  to  the  glass  slide 
and  provides  an  ohmic  electrical  contact. 

Other  research^  indicates  that  the  conductivity  of  graphite/epoxy  com- 
posite materials  is  independent  of  frequency  from  DC  to  50  MHz.  Because  the 
electrical  behavior  of  composites  is  governed  by  the  electrical  conductivity 
of  the  fibers,  electrical  examination  of  the  fibers  at  DC  using  a voltage 
source-voltmeter-ammeter  combination  or  at  400  Hz  using  a standard  curve 
tracer  provided  sufficient  information.  The  fiber  condi ctance,  G,  for  low- 
level  excitation  is  simply  the  ratio  of  the  current  I to  the  voltage,  V.  With 
this  ratio,  the  conductivity  o for  current  flow  along  the  fiber  axis  may  be 
calculated  from  a knowledge  of  fiber  length  L and  cross-sectional  area  A, 

G -I-5A 
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To  further  characterize  the  fibers  physical  1 v , both  optical  and  scanning 
electron  micrographs  of  the  individual  fibers  were  made.  Some  fibers,  as  re- 
ceived, were  partially  covered  with  a white  ash-like  substance  but  tills  was 
readily  removed  when  the  fibers  were  immersed  in  methanol. 

All  measurements  were  taken  with  the  fibers  in  room  atmosphere  with  a 
o 

temperature  of  20  C and  a relative  humidity  on  the  order  of  60%.  In  some 
situations,  it  was  useful  to  vary  the  temperature  of  the  fibers  and  tills  was 
accompl ishod  by  means  of  a resistance  heated  air  gun. 

Tile  fibers  measured  had  an  average  DC  conductivity  of  20,000  mhos/m  with 
a range  of  14,000  to  30,000  mlios/m.  All  of  the  fibers  displayed  linear  elec- 
trical responses  up  to  maximum  electric  Helds  of  approximately  4000  volts/m. 
Thirty  percent  of  the  samples  broke  at  electric  fields  below  tills  value. 

In  some  geographic  configurations  graphite/epoxy  composites  display  non- 
linear electrical  character  1st ics  at  much  lower  values  of  electric  field.  To 
further  understand  t>  » difference  between  the  bare  fibers  and  the  graphite/ 
epoxy,  a number  of  graphite  fibers  were  coated  with  epox".  These  fibers  had 
lower  threshold  power  densities.  If  the  ambient  temperatures  were  increased 
via  heat  gun,  the  epoxy  coated  fibers  would  break  at  excitation  levels  at 
which  they  were  stable  at  room  temperature.  In  one  case,  the  electric  field 
across  the  sample  was  held  at  3500  volts/m  and  the  current  remained  constant. 
Upon  application  of  the  external  heat,  the  current  increased  and  the  fiber 
failed.  The  current  increase  observed  would  seem  to  rule  out  the  force  of 
the  heat  gun's  air  stream  as  the  cause  of  fiber  breaking. 

The  detailed  electric  field  and  current  densiLy  thresholds  associated 
with  breakdown  are  shown  in  Table  1 below. 


Breakdown  Thresholds 


E(volts/m) 

minimum  3,2(10^) 

average  3.7(10^) 

maximum  4.4(10^) 

TABLE  I 


2 

J (amp/m  ) 
57.6(106) 
103.2C106) 
125.1 (106) 


These  results  Indicate  that  thermal  mechanisms  are  involved  in  the  fail- 
ure of  the  graphite  fibers.  However,  the  modeling  required  to  analyze  the 
situation  is  complex.  The  energy  generated  in  the  fiber  via  resist ive  heating 
will  be,  in  steady  state,  balanced  by  three  processes:  heat  conduction  into 
the  glass  slide,  heat  convection  into  the  air  and  radiation.  No  models  for 
conductive  heat  flow  from  a cylinder  into  a rectangular  glass  plate  could  be 
found  with  the  desired  attributes  of  accuracy  and  computational  simplicity. 

As  mentioned  above,  no  visible  color  changes  in  the  fiber  were  observed  which 
would  have  indicated  increases  in  temperature  as  large  as  lOOO^C.  For  this 
reason,  it  is  likely  that  local  temperatures  increase  at  points  along  the 
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fiber,  possibly  due  to  the  presence  of  local  constrictions  in  the  fiber  cross- 
sectional  area.  These  defects  give  rise  to  locally  enhanced  temperatures  and 
the  resultant  incandescent  failure  which  we  observed. 


Figure  4 illustrates  a simple  model  in  which  a local  defect  of  length  6 
decreases  the  area  from  A to  A'.  The  DC  current  must  be  constant  along  the 

fiber, 

J ' A*  - JA,  E'  - J'/o,  E - J/o, 

where  primed  quantities  denote  those  in  the  constricted  area. 

Define  n **  A/A',  J'  » nJ(  E'  “ HE 

The  power  density  P is  the  produce  of  current  density  and  electric 
field. 

2 

p'  * n p 

A 30£  reduction  in  area  will  result  in  a doubling  of  the  local  power 
density.  This  naturally  produces  a rise  in  fiber  temperature  with  an  accom- 
panying increase  in  the  probability  of  fiber  failure  as  a result  of  oxidation. 


Scanning  electron  microscopy  revealed  significant  surface  defects  but  it 
was  not  possible  to  determine  the  restriction  in  areas  associated  with  the 
defect  structure.  Numerical  estimates  of  the  local  heating  of  the  fibers  can- 
not be  made  without  this  information, 

The  relationship  between  tills  median ' sm  and  lightning  degradation  is 
clear  although  quite  difficult  to  model  analytically.  When  a large  current 
is  injected  into  a grapliite/epoxy  composite,  the  several  orders  of  magnitude 
difference  between  the  conductivities  of  the  graphite  and  the  epoxy  result 
in  the  total  current  being  confined  to  the  fibers.  In  the  vicinity  of  the 
attachment  point,  these  currents  will  certainly  exceed  the  threshold  currents 
for  fiber  failure.  As  fibers  break,  the  current  is  forced  into  other  fibers 
and  more  breakage  occurs.  Associated  with  this  is  epoxy  softening  and  mech- 
anical tailure. 

In  summary,  although  individual  graphite  fibers  are  difficult  to  separate 
from  tows  and  handle,  they  are  easily  characterized  electrically.  Their 
current-voltage  characteristics  arc  linear  until  the  fibers  break  and  the 
fibers  have  an  average  conductivity  of  2(10^)  mhos/meter. 


5)  BORON  FIBERS 

Boron  fibers  were  supplied  by  Avco  in  three  different  diameters:  4,  5.6 
and  8 mils  (0.01,  0.14  and  0.02  cm).  Th"se  fibers  are  produced  by  chemically 
depositing  boron  on  a heated  tungsten  wire  of  0,0018  cm.  diameter  in  an  atmos- 
phere of  ”2  + SCi^.  A reduCi. iuu  Lead  ion  occurs, 

3»2  + 2BC13  - 2b  + 6HC1 

The  boron  is  deposited  on  the  heated  tungsten  wire  and  complete  reaction 
of  the  tungsten  occurs.  The  remaining  core  consists  of  a mixture  of  WB^  and 
W B . The  core  diameter  is  unchanged  and  is  essentially  coated  with  pure 

C ->  l T 

boron. 1 


f 1 OW 


A1 1 of 

mb  or) 


the  processing  details  (exact  temperature, 
arc  proprietary  and  are  not  available- 


gas  concentrations  and 


In  order  to  characterize  the  as-grown  fibers  electrically,  it  was  neces- 
sary to  make  ohmic  electrical  contact  and  this  proved  to  be  challenging.  A 
procedure  identical  to  that  used  for  graphite  was  tried  initially.  Boron 
fibers  were  cleaned  via  ultrasonic  agitation  in  ti ichlorethylene  and  rinsed 
in  acetone  and  methanol.  The  fibers  were  dried  in  prepurified  nitrogen  and 
mounted  on  microscope  slides  using  a conductive  silver  ink.  Care  was  taken 
to  insure  the  isolation  of  the  fiber  ends  so  as  to  force  current  to  flow 
through  the  boron  sheath  radially  inward  as  shown  in  Figure  5.  The  resultant 
current-voltage  characteristics  are  sketched  in  Figure  6.  Not  only  are  they 
non-linear  but  they  display  memory. 

Tr  was  evident  that  these  results  could  not  be  used  tc  cal  ulatc  nr,  ir. 
trinsic  conductivity  of  the  boron  fibers.  Variations  of  the  cleaning  procedure 
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I-V  CHARACTERISTICS  OF  BORON  FIBERS  WITH  INDIUM  CONTACTS 
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were  tried  in  attempts  to  etch  any  contaminated  surfaces  which  could  have 
been  responsible  for  non-linear  effects.  These  treatments  included  etches 
in  both  nitric  and  sulfuric  acid.  The  rectifying  nature  of  the  silver-boron 
interface  and  the  non-linear  current-voltage  characteristics  persisted  al- 
though various  acid  concentrations  and  etching  times  were  tried. 

Contact  formation  was  next  attempted  by  "dipping"  boron  fibers  into 
molten  indium.  Relatively  pure  indium  (99.999%)  was  melted  in  a Pyrex  beaker 
over  a hot  plate  in  an  argon  atmosphere.  Cleaned  boron  fibers  were  dipped 
into  the  molten  indium  using  a glove  box.  This  procedure  naturally  resulted 
in  the  formation  of  contacts  between  the  indium  and  the  central  tungsten 
boride  core  as  well  as  the  boron  sheath.  In  situations  in  which  electrical 
isolation  was  desired  between  the  indium  contact  and  the  relatively  conduc- 
tive core,  the  ends  of  the  fibers  were  cleaved  using  a straight  edge  razor 
blade  and  the  indium  caps  removed.  Figures  7 ana  8 illustrate  this  procedure . 
It  was  during  this  cleaving  process  that  a lack  of  good  metallurgical  contact 
between  the  indium  and  the  boron  was  observable.  In  some  cases,  the  entire 
indium  contact  would  pull  away  from  the  fiber  as  the  razor  blade  was  cutting 
through. 

This  lack  of  mechanical  adhesion  was  troubling  and,  although  some  linear 
current-voltage  characteristics  were  observed  using  this  process,  instabili- 
ties and  hysteresis  persisted.  The  measured  current -voltage  curves  were  not 
necessarily  truly  indicative  of  current  through  the  fibers  only.  Rectifying 
effects  were  still  in  evidence  and  another  contacting  procedure  was  needed 
before  a conductivity  could  be  uniquely  assigned  to  the  chemically  vapor  de- 
posited boron. 

The  phase  diagrams  for  boron  and  various  metals  were  next  examined. 
Although  it  is  difficult  tc  generalize  about  the  complex  subject  of  metal- 
semiconductor  contacts,  it  is  generally  found  that  the  probability  of  fabri- 
cating an  ohmic  contact:  between  a givers  semiconductor  and  a given  metal 
increases  as  does  the  number  of  low  temperature  phases  formed  between  the 
two  constituents. 

The  various  phase  diagrams  revealed  relatively  few  candidates  among  the 
easily  deposited  metals.  Nickel  was  chosen  because  It  is  easily  electroplated 
and  forms  compounds  with  boron  at  temperatures  above  1000  C.  The  detailed 
procedures  used  for  fiber  cleaning  and  nickel  electroplating  are  given  in 
Appendix  A. 

The  nickel-boron  contacts  proved  quite  satisfactory  from  both  mechanical 
and  electrical  viewpoints.  Although  formal  tests  of  the  mechanical  adhesion 
of  the  electroplated  nickel  to  the  boron  were  not  conducted,  no  separation 
during  routine  laboratory  handling  of  the  fibers  were  observed-  In  addition, 
the  current-voltage  characteristics  of  the  boron  fibers,  with  and  without 
electrical  isolation  of  the  tungsten  boride  core,  have  been  linear. 

The  geometry  used  in  the  i-v  measurements  is  shown  in  Figure  9.  In  order 
to  relate  the  measured  conductance  to  the  boron  conductivity;  the  conductive 
core  was  isolated  from  the  nickel  plating  using  a suitable  wax  (Appendix  A). 
To  model  current  flow  in  the  fibers,  it  is  assumed  initially  that  the  core 
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conductance  is  large  compared  to  that  associated  with  the  boron  and  voltage 
drops  along  the  boride  core  are  neglected.  Further,  the  contacts  are  circu- 
larly symmetric  and  this  implies  a circularly  symmetric  electrostatic  poten- 
tial. For  the  geometry  shown  in  Figure  9,  we  write  the  static  potential  as 

<$>(r)  3 B £n(a/r) 

where  a is  the  radius  of  the  boride  core, 

R the  outer  radius  of  the  boron  sheath  and 
B is  a constant  fixed  by  the  boundary  conditions. 


CEOMETR'.  FOR  CALCULATION  OF  0„ 

O 

FIGURE  9 


By  symmetry,  a voltage  drop  of  V/2  will  appear  across  each  nickel  coated 
boron  segment  when  a voltage  of  V is  applied  across  the  entire  fiber.  Using 
this  fact,  the  constant  B in  the  above  equation  can  be  evaluated.  The  total 
potential  takes  the  form. 


Cr ) 


V £n(a/r) 
2 £n(a/R) 


From  this,  it  is  straightforward  to  calculate  the  electric  field  and  the 
current  density  and  finally'  the  total  current  I. 

„ . , tr  o V L ,1,  ,-£r.(R/a) x 

1 ~ -27IrLJ  " IKuTr)  or  0 “ (v>  <“1Tl~) 

(the  minus  sign  is  needed  on  the  left  since  positive  J implies  current  flowing  out) 
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The  first  measurements  were  made  on  samples  with  8 mils  diameter  and 
yielded  the  following  results. 


Sample  Number 

V (volts) 

I (amps) 

-3 

o (mhos/m) 

AA  6 

0.  5 

.54(10  ) 

-3 

.0837 

AA  7 

1.0 

1.05(10  ) 

, -3 

.0814 

AA  8 

1.0 

1.05(10  ) 

.0814 

(for  all  samples,  L = 1 cm.). 


In  all  cases,  the  current-voltage  curve  was  linear  for  both  polarities 
of  applied  voltage  and  displayed  no  evidence  of  any  rectifying  behavior  at 
the  nickel-boron  interface.  In  addition,  the  calculated  values  of  conduc- 
tivity are  consistent  from  f tber-to-f iber . 

Subsequent  measurement  of  a total  of  seven  samples  taken  from  the  same 
spool  of  fiber  yielded  essentially  identical  current-voltage  curves  and 
resulted  In  an  average  boron  conductivity  of 

o = 0.093  mhos/m. 

Next,  four  fibers  of  5.6  mil  diameter  were  nickel  plated  and  the  conduc- 
tivity calculated  using  the  model  developed  above.  The  results  of  this 
analysis  were 


Sample  Number 


(mhos/m) 


BB  1 


0.1 


BB  2 


0.155 


BB  3 0.097 

bb  4 0.1 


o = 0.113  mhos/m. 


Although  there  is  more  scatter  in  these  values  than  the  corresponding  8 
mil  diameter  fiber  values,  the  cluster  is  good  and  yields  a comparable  value 
for  the  average  conductivity  of  the  boron  sheath.  It  is  likely  that  thick- 
ness variations  in  the  boron  sheath  are  greater  in  5.6  mil  fibers  than  in  8.0 
mil  fibers  and  produce  larger  conductance  scatter. 

The  above  set  of  experiments  provided  initial  electrical  characteriza- 
tion of  the  "as-grown"  boron  fibers.  In  order  to  more  fully  understand  the 
mechanism  of  current  flow  in  the?T  fibers,  a series  of  electrical  measure- 
ments were  made  on  fibers  in  which  nickei  contacts  were  plated  directly  to 
the  tungsten  boride  core  as  well  as  the  boron  sheath.  The  fibers  were  5.6 
mils  in  diameter  and  typical  results  are  shown  in  the  following  table. 

24 


h«*^*tH****WM, , 





Sample 

Length 

(cm) 

Plated  length 
(each  end) 
(cm) 

Direct  Contact 
to  Core? 

DC 

Resistance 

(ohms) 

1 

9.0 

1.0 

No 

800 

2 

9.0 

1.0 

Yes 

300 

3 

16.0 

1.0 

No 

1120 

4 

16.0 

1.0 

Yes 

500 

5 

2.8 

1.0 

No 

714 

These  five  samples  are  sketched  in  Figure  10. 
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SAMPLE  GEOMETRIES 
FIGURE  10 


It  is  assumed  that  the  core  is  homogeneous  with  respect  to  fiber  length  f 

and  the  core  is  characterized  by  a resistance  per  unit  length,  r (ohms/cm).  ! 

When  this  value  is  determined,  the  electrical  conductivity  of  thS  core  will  j 

be  fixed  because  the  core  area  is  known.  The  symbol  Rg  is  used  to  denote  the  I 

resistance  of  a nickel  coated  boron  cylinder  of  lsngth°l  cm.  This  is  the  ! 

resistance  associated  with  a radial  current  from  the  nickel  to  the  conductive 
core.  In  terms  of  these  symbols,  expressions  for  the  resistances  of  each  of 
these  five  samples  may  be  written. 

The  factor  of  2 before  RR  is  a reflection  of  the  existence  of  two  radial 
paths  tliLougii  the  boron  sneatn  in  each  sample.  1 
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Sample 

1 


Expression  for  Resistance 


8r  + 2R„ 
c B 


2 9r 

c 

3 15  r + 2R 

C B 


16r 

c 


5 


1.  8r 

c 


+ 2RB 


The  algorithm  used  to  write  these  equations  should  be  clear.  The  quan- 
tity Rg  does  not  appear  in  the  expressions  for  Samples  2 and  4 because  the 
resistance  of  the  boron  sheath  is  shorted  by  the  direct  contact  between  the 
nickel  and  the  core.  Electrical  isolation  in  Samples  1,  3 and  5 was  provided 
by  masking  the  ends  of  the  boron  fibers  with  a non-conduct ive  wax  prior  to 
the  electroplating. 

Using  these  five  expressions  and  the  measured  values  of  resistance, 
values  for  r and  R can  be  obtained.  The  measured  values  of  resistance  for 
Samples  2 anS  4 lead  directly  to 

2 {9  cm.)  r = 33.3  ohms/cm 

c 

4  (16  cm)  r “31.3  ohms/cm 

c 

The  average  value  of  the  core  resistance  obtained  for  all  of  the  samples 
measured  during  the  research  period  was  32  ohms/cm.  Using  this  average  value, 
the  quantity  R may  be  calculated  using  the  experimental  results  from  Samples 
1,  3 and  5. 


1 (9  cm. ) 

3 (16  cm.) 

5 (2.8  cm.) 


800  - (8) (32) 

2 

1120  - (15) (32) 
2 

714  - (1 .8) (32) 

2 


272  ohms 

” 320  ohms 

* 328  ohms. 


There  are,  of  course,  variations  in  this  value  caused  by  differences  in 
the  "as  grown"  fibers  as  well  as  uncontrollable  variations  in  the  total  length 
of  the  plated  electrodes.  The  averaging  of  all  available  data  yields  a value 

of  300  ohms  for  R . 

'■  B 
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Thi3  result  allows  a more  accurate  calculation  of  the  boron  conductivity 
than  was  possible  in  the  last  section  in  which  the  resistance  of  the  tungsten 
boride  core  was  assumed  to  be  zero.  The  value  obtained  is 


0^  “ .25  mhos/m. 

This  value  is  two  to  three  orders  of  magnitude  greater  than  that  reported 
by  other  investigators  for  pure,  single  crystal  boron.  This  discrepancy  is 
not  alarming  in  view  of  the  large  differences  between  the  fibers  and  the 
materials  used  by  other  investigators.  In  particular,  the  boron  of  which  the 
fibers  are  made  is  polycrystalline  snd  the  impurity  content  is  simply  not 
known.  As  mentioned  above,  the  details  of  growth  are  proprietary  information. 
Because  boron  is  a semiconductor,  the  presence  c>f  relatively  small  amounts  of 
impurities  would  account  for  a large  change  in  conductivity  as  discussed  in 
Section  3 of  this  report. 

For  a fiber  of  0.01  cm.  diameter,  it  is  possible  to  define  an  effective 
conductivity  for  current  along  the  fiber.  The  fiber  conductance  G is 
written 


where  A is  the  fiber  area  and  L the  length.  But  this  conductance  is  the 
sum  of  the  conductances  of  the  core  and  the  boron  sheath. 


G - G + G 
c 

o A 
c c 

m 

L 


B 

+ 


B B 


where  the  subscript  c denotes  core  quantities  and  B implies  boron  sheath 
quantities. 

The  effective  fiber  conductivity  can  be  written  as 


o 


eff 


oc(Ac/A)  + °b<Ab//.>. 


All  three  areas  and  og  are  known.  The  core  conductivity  can  be  readily 
calculated  from  the  known  value  of  rcslstance/length,  i.e. , 

r “ l/o  A or  o ■ 1/r  A 
c c c c c c 

The  core  diameter  is  C.0018  cm,  r ■ 32  ohras/cm  and  the  core  conductivity  is 

c 

calculated  as 

o »3  (10^)  mhos/m. 
c 
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The  effective  fiber  conductivity  is  then  determined. 


o ■ 2.3  (.1.0^)  mhos/m. 
etr 


Essentially  all  of  the  current  is  confined  to  the  highly  conductive  core. 

The  above  information  concerning  the  electrical  conductivity  of  graphite 
and  boron  fibers  was  gathered  at  room  temperature,  300  K.  It  is  of  interest 
to  investigate  the  temperature  dependence  of  the  conductivity  because  this 
behavior  provides  information  as  to  the  relative  amounts  of  impurities  present 
in  a given  semiconductor  sample.  It  is  well  known  in  semiconductor  device 
physics  that  such  materials  display  a conductivity  versus  temperature  curve^ 
similar  to  that  shown  in  Figure  11.  It  is  important  to  recognize  that  the 
logarithm  of  the  conductivity  is  plotted  versus  the  reciprocal  of  the  abso- 
lute temperature  T.  This  means,  of  course,  that  temperature  increases  as  one 
moves  from  right  to  left  along  the  horizontal  axis. 

There  are  three  distinct  regions  associated  with  this  curve.  The  first, 
occurring  at  high  temperatures,  is  termed  the  intrinsic  region  and  the  conduc- 
tivity displays  an  exponential  dependence  on  the  reciprocal  of  temperature. 
Although  the  theoretical  details  are  not  of  interest  here,  the  physical  reason 
for  this  behavior  is  simply  that  the  thermally  generated  electrons  and  holes 
increase  in  number  as  the  temperature  increases  and  finally  attain  such  large 
values  that  they  swamp  out  any  effects  associated  with  electrically  active 
impurities.  All  semiconductors  will  eventually  attain  this  intrinsic  behavior 
as  their  temperature  is  increased. 

The  second  segment,  the  so-called  extrinsic  region,  is  one  in  which  the 
thermally  generated  holes  and  electrons  are  insignificant  in  number  in  com- 
parison with  those  associated  with  donors  and/or  acceptors.  In  this  region, 
changes  in  temperature  still  affect  the  thermally  generated  carriers  but, 
since  they  are  present  in  negligible  quantities,  this  has  no  effect  on  the 
sample  conductivity.  Finally,  at  very  low  temperatures,  a third  regime  be- 
comes apparent  in  which  the  conductivity  decreases  with  decreasing  temperature 
because  the  electrical  effects  of  the  donors  and  the  acceptors  are  vanishing 
as  the  electrons  and  holes  associated  with  these  impurities  "freeze"  back 
onto  their  parent  impurity  atoms.  This  process  normally  begins  at  tempera- 
tures in  the  vicinity  of  100  K and  is  not  considered  significant  in  our 
temperature  investigations. 

No  detailed  studies  of  conductivity  variations  with  temperature  were 
performed  but  a series  of  preliminary  measurements  were  made  using  0.02  cm. 
diameter  boron  fibers.  The  plated  nickel  contacts  were  isolated  from  the 
sample  ends  and,  consequently,  from  the  highly  conductive  boron  tungstide 
core. 


The  fiber  was  mounted  on  a glass  slide  which  was  placed  atop  a large 
aluminum  block  on  a hot  plate.  The  block  provided  additional  thermal  mass 
tfuu  prevented  rapid  iluuluttLiuiia  lii  Lcwpci'aLui'c  <xb  Cue  livC  plate  licdtlug 
elements  cycled  on  and  off,  A chromel-alumel  thermocouple  was  placed 
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adjacent  to  the  boron  fiber  and  the  voltage  developed  across  this  junction 
when  compared  to  a reference  junction  kept  in  an  ice  water  bath  was  measured 
on  a potentiometer. 


In  o 


1/T 


conductivity-temperature  profile 
picure  II 


This  fiber  resistance  R is  shown  as  a function  of  the  temperature 
T (°C)  in  Figure  12. 

The  sample  resistance  R may  be  written  in  terms  of  the  resistivity  P 
as 

R - YP 

where  y is  a constant  determined  by  the  details  of  the  sample  geometry. 
From  this  it  follows  that 

•inR  * £ny  + tnp 

and,  neglecting  thermal  expansion  effects, 

3-g.nR  3-Enp  dinrr 

3 (1/T)  “ 3 (1/T)  “ 3 (1/T) 
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BORON  FIBER  R-T  CURVE 
FIGURE  12 


uhere  K is  effectively  a constant,  Eg  the  forbidden  energy  gap  and  k is 
Boltzmann's  constant. 

It  follows  that 


» - Eg/ 2k  . 

3(1/1) 

The  slope  of  a plot  of  £nR  vs  1/T  is  simply  the  negative  of  the  forbidden 
energy  gap  divided  by  2k.  The  values  of  Figure  13  imply  a value  of 


Eg  - 


0.19  eV 
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The  results  reveal  that  the  boron  fibers  are  in  the  intrinsic  range  at 
room  temperature  and  are  characterized  by  a forbidden  energy  gap  of  0.19  eV. 


In  R 


This  completes  the  initial  phase  of  the  low  frequency  electrical  charac- 
terization of  graphite  and  boron  fibers.  The  methods  of  contacting  alcngwlth 
the  derived  values  of  conductivity  for  both  fiber  types  provide  a benchmark 
against  which  the  results  of  various  doping  techniques  can  be  evaluated.  In 
the  next  section,  the  details  of  the  diffusion  experiments  are  presented. 

6)  DIFFUSION  EXPERIMENTS 

The  basic  equipment  used  to  Introduce  impurities  into  both  graphite  and 
boron  fibers  was  a simple  diffusion  furnace  as  sketched  in  Figure  14.  It  is 
composed  of  an  open  ended  quartz  tube  which  is  embedded  in  thermal  insulation 
and  heated  by  adjacent  carbon  resistive  heating  elements.  The  elements  are 
excited  by  60  Hz  currents  and,  when  necessary,  the  amount  of  heating  is  con- 
trolled by  a feedback  system  which  uses  as  inputs  the  desired  tube  temperature 
and  the  actual  tube  temperature  as  measured  by  a thermocouple.  In  this  simple 
furnace,  the  thermocouple  is  slid  along  the  horizontal  axis  of  the  quartz  tube 
through  one  of  the  open  ends  until  it  is  in  the  zone  in  which  the  fibers  will 


31 


32 


be  placed.  For  our  initial  experiments,  it  was  unnecessary  to  have  precise 
temperature  control  and  the  feedback  system  was  not  used  although  the  furnace 
temperature  was  constantly  monitored  using  the  thermocouple  itself.  Tempo l- 
ature  fluctuations  were  controlled  to  +50°C  manually. 

One  end  of  the  quartz  tube  is  connected  to  a fitting  which  allows  the 
introduction  of  any  gases  desired.  In  this  work,  nitrogen  and  argon  have  both 
been  used  to  provide  inert  furnace  atmospheres.  Since  the  purpose  of  the 
experiments  is  to  examine  the  effects  of  intentionally  introduced  impurities, 
it  is  essential  that  no  extraneous  impuritie  . be  allowed  access  to  the  system. 
Both  nitrogen  and  argon  are  chemically  and  electrically  inert  and  serve  tills 
function  adequately;  although  problems  did  arise  in  the  case  of  nitrogen  as 
will  be  discussed  below,  and  the  later  experiments  used  argon  exclusively. 

The  other  end  of  the  furnace  is  fitted  with  a removable  quartz  cap  which 
provides  access  to  the  tube  for  the  insertion  of  the  fibers.  The  small  size 
of  the  graphite  fibers  caused  some  difficulty  and  quartz  slides,  grooved  to 
provide  secure  locations  for  the  fibers,  were  used  as  convenient  carriers. 

As  discussed  in  Section  3,  the  addition  of  suitable  impurities  to  graph- 
ite and  boron  fibers  is  expected  to  increase  the  intrinsic  electrical  conduc- 
tivity of  these  fibers  via  a semiconductor,  doping  mechanism.  It  is  Important, 
at  this  point,  to  examine  the  possible  atomic  species  which  could  be  used  as 
electrically  active  impurities.  Both  single  crystal  and  polycrystalline  boron 
are  semiconductors  and  are  composed  of  atoms  from  the  third  group  of  the 
Periodic  Table.  Applying  standard  semiconductor  concepts,  it  ic  anticipated 
that  impurities  from  the  fourth  column  of  the  Periodic  Table  would  act  as 
donors  while  impurities  from  the  second  column  would  act  as  acceptors.  Of 
the  impurities  in  these  two  categories,  carbon  was  selected  as  an  easily 
obtainable  species.  Silicon  and  germanium  could  also  be  used  but  they  have 
higher  atomic  masses  and  are  likely  to  diffuse  more  slowly  than  carbon. 

Graphite  describes  an  allotropic  form  of  carbon  which  occupies  the  fourth 
group  of  the  Periodic  Table  and  column  five  impurities  should  be  donors  while 
column  three  impurities  should  be  acceptors.  These  impurities  are  widely  used 
in  preparing  silicon  integrated  circuits  and  are  easily  available  in  a number 
of  useful  forms.  Boron,  phosphorous,  gallium,  indium,  and  aluminum  ore  all 
possible  candidates  but  only  boron  was  used  so  as  to  reduce  the  number  of 
experimental  variables. 

Having  chosen  the  impurity  species,  a method  was  needed  whereby  they  may 
be  introduced  into  the  fibers.  Again,  Section  3 discusses  the  process  of 
thermal  diffusion  in  general  terms  without  specifying  suitable  ways  of  intro- 
ducing impurities.  There  are  two  such  processes  in  common  use:  the  first 
involves  introducing  the  impurities  via  the  gas  flow  into  the  quartz  diffusion 
tube  while  the  second  requires  that  the  fibers  be  coated  with  a suitable  com- 
pound containing  the  impurity  species  of  interest.  In  both  techniques,  the 
goal  Is  to  insure  a high  concentration  of  the  impurity  species  at  the  surface 
of  the  fiber.  It  was  considered  important  that  the  impurities  diifuse  into 
the  fibers  with  cylindrical  symmetry  so  as  to  avoid  any  complications  intro- 
duced by  angular  variations  in  the  final  impurity  profile.  For  this  reason, 
a method  was  selected  by  which  the  fibers  could  be  uniformly  coated  using  a 
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dipping  procedure.  In  any  situation  in  which  the  impurities  were  introduced 
in  the  gas  stream,  the  surface  of  the  fiber  lying  on  the  quartz  flat  would  be 
shielded  from  the  impurity  atoms.  This  is  not  a fundamental  objection  because 
the  fibers  have  small  diameters  and  uniform  impurity  distributions  would 
undoubtedly  result  from  the  second  method  if  the  diffusion  times  were  Suffi- 
ciently long  and  the  diffusion  temperatures  were  sufficiently  high  There 
was,  as  mentioned  above,  no  literature  data  to  predict  suitable  times  and 
temperatures  and  it  was  felt  that  the  coating  procedure  would  produce  the  moat 
uniform  results  in  this  set  of  initial  investigations. 

All  fibers  were  cleaned  prior  to  coating  and  insertion  into  the  diffusion 
furnace.  The  graphite  fibers  were  cleaned  by  means  of  a three-step  solvent 
sequence;  they  were  first  immersed  in  tr ichlorethylene , rinsed  in  acetone  and 
jnethanol  and  dried  under  a heat  lamp.  The  boron  fibers  were  etched  in  hot 
sulfuric  acid,  rinsed  in  distilled  water  and  methanol  and  dried  in  a purified 
nitrogen  stream.  No  attempts  at  varying  these  cleaning  procedures  were  made. 

The  graphite  fibers  were  then  dipped  into  a commercial  borosilicate  solu- 
tion (Emulsitone  Corporation,  New  Jersey)  and  placed  under  an  infrared  heat 
lamp.  The  solution  dries  to  a glass-like  material  which  is  rich  in  boron. 

This  coating  process  is  similar  to  one  widely  used  in  the  semiconductor  industry. 

The  boron  fibers  were  coated  with  a proprietary  carbon  slurry  (CTS  Corpo- 
ration, Elkhart,  Indiana)  which  air  dries  to  an  adherent  coating.  Although 
other  carbon  sources  are  of  interest,  none  were  investigated. 

The  coated  fibers  were  placed  in  grooves  in  a quartz  slide  which  was  then 
inserted  into  the  furnace.  Of  course,  the  graphite  and  boron  fibers  were 
never  mixed  and  separate  quartz  tubes  were  used  for  each  to  avoid  cross  con- 
tamination of  impurities  which  build  up  on  the  walls  of  the  tubes. 

The  tube  endc.ap  was  then  replaced,  the  thermocouple  moved  into  position 
and  the  entire  system  flushed  with  either  nitrogen  or  argon  (each  was  used 
at  f lrst  but  argon  was  used  exclusively  in  the  later  experiments)  for  30  minutes 
before  the  resistive  heating  elements  were  excited.  The  temperature  was  then 
raised  at  a rate  limited  only  by  the  capabilities  of  the  power  supplies  and 
the  samples  were  held  at  elevated  temperatures  for  various  times. 

The  temperature  was  then  decreased  by  reducing  the  current  in  the  heating 
element  and  the  transition  to  room  temperature  normally  took  from  one  to  two 
hours.  The  inert  gas  continued  to  flow  during  this  cool  down  process  to  avoid 
the  introduction  of  impurities  at  the  surfaces  of  the  hot  fibers.  At  the  very 
least,  the  introduction  of  atmospheric  oxygen  would  result  in  some  oxidation 
of  the  fiber  surfaces  which  would  decrease  the  measured  conductivities. 

Once  the  fibers  were  removed,  the  residual  coatings  were  etched  using 
hydrofluoric  acid.  They  were  then  rinsed  in  distilled  water  and  methanol  and 
dried  in  nitrogen  and  were  then  ready  for  attachment  of  electrical  contacts 
as  described  in  Sections  4 and  5. 

The  results  of  these  diffusion  experiments  are  summarized  in  the  next 
section. 
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A brief  investigation  was  also  made  of  the  intercalation  of  graphite  fi- 
bers. This  process  is  under  investigation  at  the  University  of  Pennsylvania^ 
and  optical  data  for  intercalated  graphite  indicates  conductivities  compara- 
ble to  that  of  copper.  The  intercalation  process  involves  the  introduction 
of  large  quantities  of  super  acid  radicals  into  graphite  where  the  radicals 
act  as  electrical  acceptors  thereby  increasing  the  density  of  mobile  posi- 
tively charged  holes.  It  is  unlikely  that  this  process  would  produce  fibers 
suitable  for  advanced  composite  materials  because  of  the  corrosive  effect  of 
the  acid  upon  the  integrity  of  the  fiber-matrix  bond  which  is  critically  im- 
portant if  the  mechanical  strengths  of  composites  are  to  be  maintained. 

Following  published  procedure  , graphite  fibers  were  immersed  in  a 
mixture  of  antimony  pentaf luor ide  and  hydrofluoric  acid  for  various  times. 
Resultant  changes  in  mass  density  allowed  the  calculation  of  a stage  number. 
This  index,  in  crude  terms,  refers  to  the  percentage  of  basal  graphite  planes 
which  are  filled  with  the  super  acid  radicals.  For  example,  a stage  1 inter- 
calation implies  thac  every  basal  plane  has  accepted  as  much  acid  as  possible 
while  a stage  N intercalation  implies  that  one  out  of  every  N planes  is 
filled. 

The  first  intercalations  resulted  in  stage  numbers  in  the  range  of  10  to 
20  and,  at  such  low  acid  concentration,  no  changes  in  electrical  conductivity 
were  observed.  In  a second  set  of  experiments,  stage  5 intercalations  were 
achieved  and  observed  electrical  conductivity  increased  by  a factor  of  8 to  a 
value  of  1.6(105)  mhos/m.  Although  no  mechanical  testing  was  performed,  an 
apparently  significant  reduction  in  strength  of  the  intercalated  fibers  was 
observed  in  that  such  fibers  were  more  prone  to  break  under  normal  laboratory 
handling  than  were  untreated  fibers. 

No  significant  expenditure  of  time  was  devoted  to  the  intercalation 
process  in  view  of  the  efforts  presently  underway  at  the  University  of 
Pennsylvania  and  the  Naval  Research  Laboratory. 

7)  ELECTRICAL  BEHAVIOR  OF  DIFFUSED  FIBERS 

The  first  set  of  diffusion  experiments  used  boro3ilicate  coatings  on 
graphite  fibers  and  involved  rimes  from  one  to  twelve  hours  and  temperatures 
from  600°C  to  1200°C.  Sixteen  diffusions  were  carried  out  using  four  differ- 
ent times  (1,2,6,12  hours)  and  four  different  temperatures  (600°C,  800°C, 
1000°C  and  1200°C  ) . 

Significant  embrittlement  of  fibers  occurred  in  the  diffusion  cycles 
involving  temperatures  of  1000  and  1200°C.  These  fibers  appear  to  retain 
some  surface  residue  after  the  normal  post-diffusion  cleaning  procedures. 

In  terms  of  electrical  behavior,  the  fibers  diffused  at  temperatures  above 
800°C  display  large  resistivities  when  initially  examined.  When  a pulsed 
electric  field  of  about  three  hundred  volts  per  centimeter  was  applied,  the 
resistivity  of  the  fibers  dropped  dramatically  and  reaches  a value  of  approx- 
imately 50  percent  of  the  pre-diffusion  value.  Fibers  diffused  at  the  lower 
tcmperaturee  (600  and  800°C)  are  electrically  linear  and  display  couuucLivi- 
ties  which  are  somewhat  greater  than  the  starting  fibers  themselves.  This 
increase  in  conductivity  was  as  large  as  a factor  of  five  in  a few  fibers 
and,  in  others,  as  small  as  ten  percent. 
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Nitrogen  was  used  as  the  diffusion  furnace  gas  in  this  set  of  experiments 
and  another  set  was  carried  out  over  the  same  range  of  temperatures  and  times 
using  argon.  A borosilicate  compound  was  again  used  as  a source  of  boron  as 
described  previously.  The  results  were  essentially  equivalent  to  those  ob- 
served with  nitrogen.  Significant  fiber  embrittlement  was  observed  at  the 
highest  diffusion  temperatures,  conductivity  enhancement  was  clear  in  many 
samples,  and  some  residual  coatings  remained  after  chemical  etching. 

A literature  survey  was  conducted  in  an  attempt  to  identify  the  compound 
which  might  be  forming  on  the  surfaces  of  the  fibers.  BN^  is  a possible 
candidate  because  it  is  insoluble  in  acids-  To  check  this  hypothesis,  graph- 
ite fibers  were  Immersed  in  fused  NaCl  held  at  a temperature  of  801°C.  The 
time  of  immersion  was  10-30  seconds.  The  fibers  were  electrically  linear 
after  this  treatment  and  displayed  no  evidence  of  a surface  residue.  It  was 
felt  that  the  use  of  borosilicate  as  a source  of  boron  impurity  atoms, although 
successful  in  improving  conductivity,  leads  to  the  formation  of  chemically 
resistant  coatings  which  are  removable  in  fused  alkalis.  The  experimental 
difficulties  involved  in  this  procedure  induced  a search  for  alternate  sources 
of  boron  for  che  diffusion  experiments. 

A boron  nitride  (BN)  source  was  purchased  from  Materials  Research  Corpo- 
ration and  placed  into  the  diffusion  tube  upstream  (in  terms  of  a gas  flow) 
from  the  quartz  slide  carrying  the  graphite  fibers.  The  bcron  atoms  released 
from  the  heated  boron  nitride  are  carried  along  in  the  gas  stream  and  provide 
a boron  rich  layer  at  the  surface  of  the  fibers.  A set  of  experiments  using 
this  source  was  carried  out  at  temperatures  of  600  and  700°C  for  periods  of 
three  and  six  hours.  In  all  cases,  the  increases  in  graphite  conductivity 
were  between  30  and  50  percent. 

Another  sec  of  experiments  was  carried  cut  for  temperatures  up  to  1200°C 
for  times  ae  long  as  20  hours.  In  no  case  were  conductivity  increases  greater 
than  the  factor  5 observed. 

This  was  surprising  because  it  had  been  anticipated  that  increased  dif- 
fusion times  would  lead  to  increased  values  of  conductivity.  There  is  some 
experimental  evidence  that  the  impurity  species  will  occupy  interstitial  loca- 
tions to  a high  degree  if  the  diffusion  temperatures  are  below  2000°C.  It  is 
likely  that  the  factor  of  five  improvement  in  conductivity  which  has  been  ob- 
served is  the  natural  limit  for  impurities  diffusing  at  relatively  low  tem- 
peratures. The  disorder  in  the  graphite  fibers  may  be  such  that  substitu- 
tional positions  can  be  occupied  by  impurity  atoms  in  large  numbers  only  if 
the  diffusion  temperatures  are  high.  The  equipment  precluded  testing  this 
idea  because  the  quartz  used  as  a furnace  tube  placed  an  upper  bound  of  1300°C 
on  the  operating  temperature. 

Figure  15  summarizes  the  conductivity  enhancements  observed  in  this  first 
set  of  experiments. 

To  obtain  higher  diffusion  temperatures,  it  was  necessary  to  design  and 
construct  a small  induction  heated  furnace.  This  furnace  was  capable  of  ob- 
taining temperatures  of  2aU0°C  and  was  made  operational  in  the  last  month  of 
the  contract  period.  It  produced  encouraging  results  in  the  first  set  of 
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experiments-  These  are  summarized  in  Figure  16. 
ities  of  1(10^)  mhos/m  were  measured. 


For  some  fibers,  conductiv- 
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CONDUCTIVITY  ENHANCEMENT  IN  GRAPHITE 
FIGURE  15 


In  view  of  the  extended  difficulties  in  making  ohmic  electrical  contacts 
to  the  boron  fibers,  only  preliminary  diffusion  experiments  could  be  carried 
out.  In  all  cases,  the  quartz  diffusion  furnace  with  an  argon  atmosphere  was 
used.  The  first  experiments  were  conducted  at  diffusion  temperatures  of  1000 
and  1200°C  for  times  of  four,  twelve,  and  twenty  hours.  These  results  are 
summarized  in  Figure  17. 
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8)  SUMMARY  OF  RESULTS  AND  CONCLUSIONS 


The  activities  during  this  contract  period  have  led  to  two  significant 
conclusions.  The  conductivities  of  graphite/epoxy  and  boron/epoxy  advanced 
composite  materials  depend  primarily  upon  the  conductivities  of  the  rein- 
forcing fibers  and  the  degree  of  fiber  contact.  From  this,  it  follows  that 
the  overall  composite  conductivities  can  be  enhanced  in  a mechanically  accept- 
able fashion  only  by  increasing  the  fiber  conductivities  themselves.  Second, 
the  additions  of  boron  to  graphite  fibers  and  carbon  to  boron  fibers  result 
in  significant  increases  in  fiber  conductivities. 

The  effort  begun  this  year  will  continue  under  the  sponsorship  of  the 
Air  Force  Office  of  Scientific  Research.  The  goals  of  this  work  are  to 
develop  better  statistical  models  of  current  flow  in  advanced  composite 
materials  and  to  more  thoroughly  characterize  the  diffusion  process  itself. 

In  particular,  during  the  period  covered  by  this  report,  there  were  no  tools 
available  to  determine  the  precise  distributions  of  impurities  present  in 
the  fibers  before  and  after  the  diffusion  cycles.  One  of  the  goals  of  next 
year's  program  is  to  develop  an  inhouse  capacity  to  grow  both  graphite  and 
boron  fibers  so  that  the  background  impurity  content  can  be  minimized.  In 
addition,  the  facilities  ~f  the  National  Research  and  Resource  Facility  for 
Submicron  Structures  at  Cornell  University  will  be  used  to  measure  the  im- 
purity profiles  in  fibers.  This  will  be  done  using  a Secondary  Ion  Mass 
Spectrometer  which  uses  an  ion  beam  to  sputter  material  from  a sample  and 
then  spcctrographically  determines  the  atomic  content  of  the  sputtered 
material.  This  tool  will  determine  the  material  composition  of  fibers  as 
a function  of  radial  depth.  In  this  way,  precise  values  for  the  diffusion 
coefficients  of  impurities  in  graphite  and  boron  will  be  determined. 
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APPENDIX  A 


A-l  INTRODUCTION 

This  appendix  briefly  describes  the  technique  of  nickel  plating  as  it 
was  used  to  electrically  contact  boron  fibers. 

Ir  Section  I,  the  nickel  plating  bath  used  in  this  particular  case  is 
described.  Section  II  summarizes  the  cleaning  and  plating  procedures  and, 
in  the  conclusion,  the  success  of  this  technique  is  described. 

A- 2 THE  NICKEL  PLATING  BATH 


The  nickel  plating  solution  used  in  this  particular  experiment  is  a 
Watts  Bath.  The  main  constituents  of  this  bath  are  nickel  chloride,  nickel 
sulfate  and  boric  acid.  Nickel  sulfate  is  the  principal  source  of  nickel 
ion  in  the  Watts  Bath.  The  concentration  of  the  nickel  sulfate  determines 
the  limiting  current  density  and  the  plating  rate.  Consequently,  this 
determines  how  much  current  can  flow  through  the  cathode  (i.e. , the  boron 
fiber).  Nickel  chloride  improves  anode  (i.e.,  nickel  metal)  corrosion  and 
conductivity.  Increasing  conductivity  is  of  practical  importance  in  order 
to  reduce  the  voltage  necessary  to  produce  the  desired  current  density. 

Boric  acid  is  the  solution  buffer  and  helps  to  produce  smoother  and  more 
ductile  deposits  of  nickel. 

Sound  deposition  of  nickel  depends  on  other  parameters  besides  the 
chemicals  in  the  solution.  Temperature  agitation,  controlling  of  the  pH 
level,  anti-pitting  reagents,  optimal  anode-cathode  current  densities  and 
low  contamination  are  all  of  importance. 

The  bath  that  is  used  in  our  particular  case  has  the  following  composi- 
tion and  operation  conditions: 


Constituent 

Quantity 

Nickel  Sulfate 

300  g/1 

(NiS04.6H20) 

Nickel  Chloride 

60  g/1 

(NiCl2.6H20) 

Boric  Acid 

42  g/1 

(«3bo3) 

4%  Sodium  Hydroxide  Solution 

42  g/1 

(NaOH) 

30%  Hydrogen  Peroxide  2.5  ml/1  diluted  in  200  ml  of  Ho0 

' " 4- 

Activated  Charcoal  2.4  g/1 
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Operating  Conditions 
Temperature 

PH 

Current  density  of  anode 
Current  density  of  cathode 


Range 

49°  to  66°C 
3.5  to  5,0 


Preferred 

60°C 

4.5 

2 

100  amp/m 

2 

530  amp/m 


The  following  steps  are  typical  in  the  preparation  of  500  ml  of  solution: 


Day  1:  1,  Fill  a clean  beaker  with  500  ml  of  distilled  water. 

Heat  it  up  to  66°C. 

2.  Add  I50g  nickel  sulfate  and  30g  of  nickel  chloride. 

3.  Add  4%  sodium  hydroxide  until  the  pH  is  approximately 
5.2;  agitate  vigorously. 


4.  Add  30%  hydrogen  peroxide  - 1.25  ml  diluted  in  100  ml 
of  H20.  Stir. 

5.  Add  1.2g  of  activated  carbon,  stir  for  2 hours  and 
allow  solution  to  settle  overnight. 


Day  2:  6.  Filter  all  the  carbon  out. 

7.  Add  21g  of  boric  acid  (stirring  continuously) 

8.  Adjust  pH  to  4 with  dilute  CP  graded  sulfuric  acid. 


Finally,  for  the  first  plating  job  using  this  solution,  it  is  better  to 
use  a dummy  cathode  (e.g.  copper  wire).  Nothing  has  been  said  so  far  about 
the  anode  and  specific  characteristics  are  given  in  Section  II.  For  our  case, 
a 100%  pure  nickel  anode  is  not  necessary  which  makes  the  technique  econom- 
ically attractive. 

A- 3 CLEANING  AND  PLATING  PROCEDURES 


In  crder  to  attain  consistent  results  for  every  plating  of  boron  fibers, 
cleaning  of  the  fiber  is  of  vital  importance.  Considering  the  small  diameter 
of  these  fibers  (e.g.,  4 to  8 mils)  any  particle  of  foreign  material  at  the 
surface  of  the  fibar  will  be  covered  with  nickel  and  give  v-i  measurements 
inconsistent  with  similar  samples  of  boron  fiber.  In  order  to  avoid  this  the 
following  procedure  is  u3ed: 

(1)  Vapor  degrease  the  sample  with  trichlorethylene  for  at  least 
5 minutes. 

(2)  Run  the  fiber  through  acetone,  methanol  and  distilled  water. 
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(3)  Drop  fiber  In  sulfuric  acid  for  3 minutes. 

(4)  Run  through  distilled  water. 

(5)  Blow  dry  with  prepurified  nitrogen.  (Do  not  use  any 
paper  cleaning). 

Before  plating,  the  current  densities  were  calculated  using  the  follow- 
ing relationship: 

•*ASA  " ^CSC  " *p 

where  J ■ current  density  of  the  anode 
A 

S.  - surface  area  of  the  anode 
A 

J ■ current  density  of  cathode 
a 

S - surface  area  of  cathode 
B 

Ip  ■ plating  current 
2 

Assuming  j « .01A/cm  (preferred  value  for  best  results)  and  i » 

. 054A/cm^  the  relation  becomes  SA  ■ 5.38  S . This  immediately  determines  the 
ideal  minimum  size  of  the  nickel  metal  ano&e. 

For  our  case,  in  the  range  of  different  diameters  of  fibers  and  up  to 
20  cm  long  samples,  the  anode  is  immersed  3-4  cm.  and  the  plating  current  is 
set  to  6 ma. 

By  setting  the  temperature  at  60°  + 2°C,  the  ph  at  4.5,  T.  ■ 6 ma  and 
lightly  agitating,  best  results  are  attained.  P 


A-4  CONCLUSIONS 


The  technique  thus  far  described  has  proven  to  be  a reliable  method  of 
contacting  boron  fibers.  It  is  also  an  easy  and  economic  way  of  depositing 
uniform  sheets  of  nickel  on  the  boron  fibers.  Further,  controlled  thickness 
and  amount  of  plating  area  make  it  convenient  to  measure  v-i  characteristics 
and  verify  current  flowing  predictions  on  the  basis  of  the  sample's  geometry 
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